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Abstract 

We investigate an effective model for organometallic complexes (with potential uses in optoelectronic devices) 
via both exact diagonalisation and the configuration interaction singles (CIS) approximation. This model 
captures a number of important features of organometallic complexes, notably the sensitivity of the radiative 
decay rate to small chemical changes. We find that for large parameter ranges the CIS approximation 
accurately reproduces the low energy excitations and hence the photophysical properties of the exact solution. 
This suggests that electronic correlations do not play an important role in these complexes. This explains 
why time-dependent density functional theory works surprisingly well in these complexes. 

Keywords: effective Hamiltonian, Hartree-Fock, mean field theory, configuration interaction singles 
(CIS), organometallic 



1. Introduction 

Organometallic complexes have significant po- 
tential for use as the optically active materials in 
organic photovoltaic (OPV) devices [H-il, organic 
light emitting diodes (OLEDs) 043' organic 
light emitting field effect transistors (OLEFETs) 
U However, progress in this field has been 
hampered by difficulties in identifying useful design 
rules for new materials. The optoelectronic func- 
tionality of such complexes depends on the proper- 
ties of their lowest few excited states. Therefore, in 
order to aid with the design of the future genera- 
tion of organometallic complexes, one might think 
that models are needed to accurately describe ex- 
cited states on a case-by-case basis. However, an 
alternative approach is to introduce models that 
can describe whole classes of materials and thereby 
understand and predict trends across a range of po- 
tential complexes [loj . 

The transition metal cores; typically a iridium, 
ruthenium, platinum or palladium atom; of the 
organometallic complexes discussed in this letter 
are associated with strong electronic correlations. 
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The electrons in these materials are strongly corre- 
lated because of their tight confinement in d and / 
atomic orbitals, which increases the relative mag- 
nitude of electronic interactions ll|, |l2| and can 
have significant effects on the properties of these 
materials 13-15J. A powerful approach to mod- 
elling strongly correlated electrons is to introduce 
effective Hamiltonians via a physically-motivated 
reduction in the many-body basis. Effective Hamil- 
tonians have played a crucial role in explaining 
the str ong correlation effects in both inorganic 
H, E [lil and organic materials (iBUlTj and have 
provided new i nsig hts in to some organometallic 
compounds [13, [3 [iM [3l ■ Thus far such models 
have not been widely used to model the optical exci- 
tations of organometallic complexes with potential 
for OPV, OLED or OLEFET apphcations [13, 
22|. To date the modelling of these complexes 
has mostly been based on first principles quantum 
chemistry calculations. Effective Hamiltonian and 
first principles approaches each have their own set 
of advantages and drawbacks and combin ing the 
two approaches is often extremely powerful Ifl, 23 1. 

Because of the relatively large size of the com- 
plexes of interest for OPV, OLED or OLEFET 
applications most first principles studies are lim- 
ited to density functional theory (DFT) and time- 
dependent DFT (TDDFT). One puzzle is that 
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TD(DFT) approaches have been shown to give 
reasonably accurate predi ctions for a number of 
organometalic compounds 



24- 33t , which one would 



not expect if their are strong electronic correlations 
3J,|35[ . The complementary approach of construct- 
ing and solving effective model Hamiltonians can 
address some of the shortcomings of TDDFT. One 
of the advantages of effective Hamiltonians is that 
they can be used test the validity of approxima- 
tions and investigate the parameter regimes where 
they break down. By reducing the model to the 
few degrees of freedom needed to describe low en- 
ergy processes, one can include a more detailed de- 
scription of electronic correlations, at the cost of 
some chemical specificity. On reducing the descrip- 
tion of the system to only the effective low energy 
interactions, one often finds that the same model 
is applicable to many seemingly disparate systems. 
For example, the Anderson single impurity model 
has been applied to such varied systems as mag- 
netic impurities in metals, quantum dots in semi- 
conductor heterostructures, carbon nanotubes, and 
single molecule transistors T^, the Heisenberg 
model has found applications from magnetic sys- 
tems to quantum computers to valence bond theo- 
ries of chemistry and the Schrieffer-Su-Heeger 
model has been used to understand the properties 
of a variety of polymers 36l 37[. 

Here we investigate an effective Hamiltonian for 
heteroleptic organometallic complexes via both ex- 
act diagonalisation of the Hamiltonian (full con- 
figuration interaction) and the configuration in- 
teraction singles (CIS) approximation based on 
in a Hartree-Fock reference ground state. We 
show that, for most of the parameter space of the 
model, CIS is an extremely accurate approxima- 
tion. We give criteria for when the CIS approxi- 
mation breaks down and show that CIS can accu- 
rately predict a number of trends observed in both 
the exact solution of the model and in experiments 
on organometallic complexes, such as large changes 
in the radiative rate from small changes in model 
parameters, corresponding to subtle chemical sub- 
stitutions on the organic ligand. 

2. Methods 

We study a model Hamiltonian for complexes 
where one metal orbital and one pair of frontier lig- 
and orbitals dominate the low energy physics [22| 
(this is observed in a wide variety of heterole" 
tic complexes for example Ru(NH3)2Cl2(bqd: 



Ru(dcbpy)(bpy)2 128|, or Pt(cnpmic) [21|, further, 
it has been argued that solvent effects or excited 
state geometry relaxation can cause one ligand to 
dominate the physics of homoleptic complexes such 
as Ir(ppy)3 [11]): 
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where fin = Y^a^lH^, flL = Y^a^a^a, 

J2a 13 ^a'^al3^l3' ^'^^ annihilates (creates) an elec- 
tron in the 'renormalized HOMO level of the lig- 
and', Lct'' annihilates (creates) an electron in the 
'renormalized LUMO level of the ligand', m'!^^ an- 
nihilates (creates) an electron in the 'renormalized 
metal orbital', e is the difference in the renormal- 
ized energies of the metal orbital and the ligand 
HOMO, A is the renormalized HOMO-LUMO gap, 
(t^) is the effective hopping amplitude between 
the metal orbital and the Hgand HOMO (LUMO), 
J is the effective exchange interaction between the 
ligand HOMO and the ligand LUMO, is the ef- 
fective Coulomb repulsion between two electrons 
in orbital i {i £ {H, L, M}), and Vij is the effec- 
tive Coulomb repulsion between an electron in or- 
bital i and another in orbital j. The single elec- 
tron orbitals are chosen to be orthonormal, i.e., 
{H.,M,} = {Ha.L,} = {Ma,L,} = 0. For the 
discussion of the solution of this model, below, it 
will be helpful to define the state |uac), the 'vac- 
uum' state with {fiu) — {nu) = {fih) = 0. 

It is important to note that renormalized fron- 
tier (HOMO, LUMO and metal) orbitals may be 
quite different from the bare HOMO, LUMO and 
metal orbitals. The state H'l\vac) (ZJ-jwac)) will, 
in general, be significantly different from the high- 
est occupied (lowest unoccupied) molecular orbital 
of the ligand found via a Hartree-Fock calculation 
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on an isolated ligand, and similarly the 'metal or- 
bital' Enight have electron density on the ligands 
lol . 22, 33- This means that the second quantized 
notation we adopt in this letter is significantly more 
convenient than wavefunction based notations. 

The set of interactions retained in Hamiltonian 
([T]) is equivalent to an INDO (incomplete neglect 
of differential overlap) Hamiltonian 34^, in that we 
have included all direct Coloumb interactions, but 
only the local exchange interactions. Like the one 
electron orbitals, the Hamiltonian parameters are 
renormalized. In principle one could derive the pa- 
rameters by integrating out high energy degrees of 
freedom from first principles calculations 3^4^. 
However, this is a highly computationally demand- 
ing task therefore, a more fruitful approach is to 
make use of semi-emprical parameterisations 



To find the Hartree-Fock ground state of this 
Hamiltonian one constructs a variational single 
Slater determinant state and minimizes the energy 
of that state To construct the single Slater 

determinant, we first chose, without loss of gener- 
ality, an arbitrary set of three mutually orthogonal 
one-hole operators, corresponding to linear super- 
positions of the three orbitals in this model 

do- = sin (f) sin 6 Ha + cos 

+ COS0Lcr, (2) 

ha = sin (3 ^cos (f) [ sin 9Ha + cos 9M, 

— sin (pLa 

+ cos l3 (^cos 9 Ha -sine Ma^, (3) 
Ca = COS /3 ^ COS (/) [ sin 9 Ha + cos 9Ma] 

— sin (j)La 

- sin (3 cos 9Ha - sin 9Ma^ . (4) 

These states are orthonormal for any choice of the 
angles 9,(j) and /3, e.g. {aa,ba} = {ba,Ca} = 
{Ca,aa} — (since, by definition, iJ, M and L 
act on orthogonal spaces). 

We construct a closed shell four electron single- 
determinant ground state from one of these arbi- 
trary states, writing it as two holes in the six elec- 
tron state 
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a^a-f (^lIlImImIhIhI^ \vac) 
a4,atl*6) (5) 



cos^ 0|O) -I- sin^ ^sin^ (j)L\LlH^H^\0) 
+ cos^ 9 sin^ 0i|i|M^A/ft|O) 
cos 9 cos 4> sin t/i 

sin 9 cos (f> sin 
1 

72 

cos 9 sin 9 sin (h 

x^l\lI{H^M^-H^M^)\Q) (6) 

cos^ 0|O) + sin^ 6'sin2 ^l^iC^) 
-hcos^ 6* sin^c^l ^MLCT^) 
+V2 cos 9 cos <t> sin (t>\^MLCT^) 



x—{L\H^ + LlH^)\0) 



-V2sm9 



sm cos sm c 



- V2 cos 61 sin 9 sin^ (/)| ^ MH^) , 



(7) 



where |0) = M^U/fj^ijt^t|„ac) = i^^i^tl^e 
I $6) is the state with 6 electrons, i.e. every or- 
bital doubly occupied. It is worth briefly discussing 
the basis in which we have expressed the varia- 
tional state in Eqs. ([6]) and ([7]), as it provides a 
natural basis in which to study this system and 
we will refer to these states frequently. This ba- 
sis of states comprises of the singlets defined by 
Eqs. and ([7]) and their triplet counterparts has 
{fin) T {tlm) and (fji) as quantum numbers. These 
states are the eigenstates of the Hamiltonian, Eq. 
(IT|), for = = 0. The superscript prefix in the 
state label denotes spin degeneracy, and the super- 
script suffix denotes LUMO occupancy {ul). The 
reference state, |0), has no electrons in the ligand 
LUMO, and two each in the HOMO and metal or- 
bitals. \^LC^) and \^LC^) are, respectively, singlet 
and triplet ligand centred excitations, where one 
of the HOMO electrons of the reference state has 
been excited to the LUMO. In the metal-to- ligand 
charge transfer states, \^MLCT^) and \^MLCT'^), 
an electron has been moved from the metal orbital 
to the LUMO. \^LC^) and \^MLCT'^) are the states 
where two such ligand centred or MLCT excita- 
tions, respectively, have occurred, while |^A/_ff^) 
and \^MH'^) have one excitation of each kind. 

Eqs. (HI and © give a natural interpretation 
of the angles 9 and </>. (j) determines the degree of 
ligand LUMO character in the ground state. For 
(j) — the ground state is pure reference state with 
{ul) — 0, whereas for </> = 7r/2 the ground state 



i^L-|-|$fi) and 
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has (ni) = 2. The angle determines the degree of 
MLCT or LC character in the ground state, with 
6* = meaning purely MLCT {{uh) = 2) and = 
tt/2 purely LC {{um) = 2). From Eqs. and g]) 
we can see that /3 determines the degree of ligand 
LUMO versus HOMO-metal character in the b- and 
c-hole states. 

The energy of the Hartree-Fock state is given by 

= ($6l44^^-L"tl*6), (8) 

and must be minimized with respect to 6 and </>. 
This minimum is the Hartree-Fock ground state en- 
ergy E^'', defined as E^^ = E^^iOa, (po)- 

For the CIS calculation we need to find the singly 
excited states. The singly excited states are con- 
structed by replacing one of the a-holes with a b- or 
c-hole. The Sz ~ ±1 states are of the form 

haUsD=b,aUiat\^6) (9) 

or 

c,at|5o^^) ^c^aU^atl^e) (10) 
while the Sz = states are of the form 

-i= (b^al ± 6^4) a4.atl*6> 

= ^ (^iatTMi) 1*6), (11) 

or ^ 

— (c_^4±ct4) a^atl^e) 

= ^(c;atTcta;)|$6), (12) 

which ensures that they are spin eigenstates. The 
Sz — states defined above are a combination of 
two Slater determinants, while the Sz — O- states 
are single Slater determinants. To perform the CIS 
calculation one computes the Hamiltonian matrix 
in this basis of singly excited states and solves it 
self-consistently [44,^]. In this case, the remaining 
variational angle f3 allows us to directly minimize 
one of the singly excited states energies, equivalent 
to diagonalizing the CIS Hamiltonian matrix. Here 
we choose the singlet excited state defined in Eq. 
(ITlT) as the lowest excited state. 

We perform the minimization of the lowest sin- 
glet and triplet states with respect to (3 separately 
for the singlet and triplet sectors (leading to a dif- 
ferent value of /3 for each spin sector, which we de- 
note f3s and /3t respectively). Thus all three free 



parameters in the definition of the one-hole basis 
states are fixed for both the singlet or triplet sectors 
([6*0, 00, /3s] and [0o,0o,/3t], respectively). Hence, 
the energy of every excited state can now be com- 
puted. 

3. Results 

We now explicitly compute the CIS energies 
and compare them to the exact eigenvalues of the 
model. To do this we use a standard set of param- 
eter values around which we will vary one param- 
eter at a time to investigate its effects. A typical 
set of parameters for organometallic complexes is: 
J = 1 eV, A = 3 eV, £ = 0.25 eV, t" = 0.1 eV, 
= 0.1 eV, Um = Uh = Ul = U = 3 eV and 
Vhl = Vhm = Vlm = V = 3 eV UM- 

We find that the CIS excited states energies are 
typically extremely close to the exact energies. Fig. 
[1] shows that in the ground state and singly excited 
states the error in the energy is well below 1%. It is 
also interesting to compare both the exact and CIS 
results to the (exact) solutions for = = 0. 
In this limit the solutions are eigenstates of hn, 
fiM and fiL and thus are given by the basis defined 
by and below Eqs. © and ©. When t" Q ^ 
these basis states mix. Significant mixing only 
occurs when the energy difference between the — 

— eigenstates is the same order as the matrix 
elements (proportional to or t^) coupling them. 
Since and are expected to be quite small (0.1 
eV in the parameterization used here), for most of 
parameter space the exact eigenstates are close to 
being pure = ^ eigenstates. 

Fig. [IJalso shows that around the avoided cross- 
ings, the error in the higher excited states ener- 
gies increases (even if the crossing does not involve 
any of the high lying excited states, such as near 
e = —0.75 eV). Recall that 9 and 4> are determined 
by optimizing the ground states energy, and (3 is de- 
termined by optimizing the first excited states en- 
ergy. Thus, once the energies of these states have 
been optimized, there are no more free parameters 
to optimize the higher excited states. Neverthe- 
less, away from avoided crossings the CIS calcula- 
tion does a remarkably good job of predicting all 
the the excitation energies, even those of the dou- 
ble excitations. 

Fig. [5] shows that this is also true of the triplet 
excited states. Once the ground state energy is op- 
timized there is just one minimization parameter 
left in the triplet subspace, /3t , and yet all three 
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triplet states are within 0.1% of the exact energy 
over the range plotted here. 

It is important to stress that having a good 
variational estimate of a states energy is a neces- 
sary but not sufficient condition for having a good 
variational wavefunction. Accurate wavefunctions 
are necessary, however, to make accurate predic- 
tions of the physical and chemical properties of 
the complexes under consideration. Fig. [3] shows 
the error in the singlet wavefunctions, defined as 
1 - I (S-f I Sf I . It is clear from this figure that 
around the avoided crossings at e = —0.75 eV, 
eV, and 1.1 eV the higher excited states are not cor- 
rectly reproduced by the mean field CIS procedure. 
Fig. |3] also makes it clear that the singly excited 
states are much closer to the exact solutions than 
the doubly excited states. 

Further understanding of the successes and fail- 
ures of the CIS approximation can be gained by 
examining the wavefunctions of the excited states. 
For example. Fig. 2] shows the dominant compo- 
nents of the exact and CIS solutions for the third 
excited singlet state, S3. Note that, except near 
the avoided crossings at e = eV and 1.1 eV, 
one — — eigenstate dominates the ex- 
act eigenstate. The singlet eigenstates are nearly 
pure = = eigenstates because the hop- 
ping terms and mixing these states are ex- 
pected to be small compared to the other Hamilto- 
nian parameters, which determine the energy gaps 
between these states. The = ^ eigen- 
states are spin symmetrised pairs of Slater determi- 
nants; these states have the minimum static correla- 
tion. Near avoided crossings, when these states are 
strongly mixed in the exact eigenstates, the exact 
states have strong static correlations. The CIS solu- 
tion does a poor job of reproducing the exact state 
near the avoided crossings at e = —0.75 eV, eV 
and 1.1 eV, as seen previously in Fig. [31 In the CIS 
solution the avoided crossings are sharp crossovers 
rather than the smooth transitions seen in the exact 
solution. 

When constructing the variational ground state 
Eq. ([5]) it was assumed that the exact ground state 
is closed shell. We see in Fig. [S]that if the Coloumb 
energy Um becomes large then this assumption is 
not valid. Above U]\i = 6 eV the exact ground state 
is predominantly the open shell \^MLCT^) state. 
Since there is no way to write an open shell state in 
terms of two a-holes the error in the wavefunction 
increases quickly, and is over 25% for Um — 6.5 eV. 
Thus, this model has a clear dividing line between 



a weakly correlated regime, where the ground state 
is closed shell, and a strongly correlated regime, 
where the ground state is open sheU. The CIS ap- 
proximation fails dramatically when in the strongly 
correlated regime {Um ^ 5.5 for the parameters in 
Fig. [5]). However, one does not expect the ground 
state to be open shell in organometallic complexes 
and thus the CIS approximation should give reliable 
results for these materials. 

Now that we have have established the regime 
of validity for the CIS approximation, we can ap- 
ply it to calculate other properties of the excited 
states with the goal of understanding the radiative 
properties of organometallic complexes. Of partic- 
ular interest is the phosphorescent lifetime of vari- 
ous organometallic complexes. Radiative emission 
from triplet states (phosphorescent decay) is made 
possible by including spin-orbit coupling. Further, 
it has been found that one must include spin-orbit 
coupling in TDDFT calculations to correctly de- 
scribe the spectra of organometallic complexes 24 1 . 
Therefore, in order to study the radiative proper- 
ties predicted by the current model we add spin- 
orbit coupling perturbatively to the solutions of the 
model Hamiltonian. A triplet state \Tm) has a tran- 
sition dipole moment, Mt^, to the ground state, 
1 5*0), that is given (to first order in the spin orbit 
coupling Hamiltonian, Hso) by 



{Tm\Hso\So) 



Es 



E. 



So 



(13) 

where Et^ is the energy of the triplet state, Es„ 
is the energy of the singlet state |S'„) with transi- 
tion dipole moment to the ground state Mg^, the 
sum runs over all excited singlet states, and the fi- 
nal term is the direct excited state to ground state 
transition (icj . 

Eq. [T3] is only valid if the effect of the spin-orbit 
coupling is perturbative, i.e., for {Tm\Hso\Sn) ^ 
Et„^ — Es^ ■ Thus we only expect the perturbative 
approach to be applicable if the all singlet excited 
states are sufficiently different in energy from the 
triplet state in question. The only singlet close in 
energy to the first triplet is the first excited sin- 
glet; all the other singlets are much further away 
in energy. In the parameter regime where the per- 
turbative method is valid, the first excited singlet 
state is dominated by the \^MLCT^) state (Fig. [T]) 
and the first triplet is dominated by \'^LC^) (Fig. 
12). This allows us to make the approximation that 
\^ M LCT^) is the only singlet state that will signif- 
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icantly couple to Ti. Further, we neglect spin-orbit 
coupling between states unless both have a large 
metal character as the magnitude of the spin-orbit 
coupling interaction scales as Z'^, where Z is the 
atomic number. Thus, Eq. (|13p simplifies to 



M 



(TilHlLCT^) 
{^MLCT^\Hso\^MLCT^). (14) 



It is helpful to write the singlet eigenstates in the 
form 



^ |o) + 4y^^^VMLCTi) 



(15) 



thus 



Ml = e{So\f\Si 



:<4o('/'l^l<^'), (16) 



where the sums over (j) and 0' run over all of the 
states. In order to evaluate the transition dipole 
moment between Si and So one then simply notes 



that 



^MLCT^) ~ 0, due to the spatial sepa- 



ration between the metal and ligand orbitals, simi- 
larly other matrix elements involving metal to lig- 
and charge transfer will be negligible. For the typ- 
ical parameters in the regime where the pertur- 
bation theory converges \Si) ~ \^MLCT^) and 
IS'o) ~ |0) (see Fig. [T]). Furthermore c%^c%^ ~ 
—Cgf^^^ Cg^^^*^"^ , which can be seen numerically 
or derived perturbatively from the — — 
limit. Hence 



Ml 



.0 „o ^^-.^ 



where Sf 
Therefore, 



Ml 



^MLCT^\f\^MLCT^) 



(17) 

if|0). 



P^O 



{Til^MLCT^) 



Eti — EiMLCT^ 

X Sr^MLCT^lHso T A/LCT^) . (18) 

From the transition dipole moment one can find 
the radiative decay rate of the triplet via the Ein- 
stein A coefficient. 



-^^\M° 



(19) 



where Eq is the permittivity of free space and huj 
is the triplet-ground state energy gap f??]. The 
quadratic dependence on Ml amplifies the effects 
of small changes in Hamiltonian parameters on the 



radiative decay rate. For concreteness, we use a 
physically reasonable value for the strength of the 
spin orbit coupling, M LCT^\Hso\^ M LCT^) = 
100 cm~^ [i^ and choose (5r such that we reproduce 
a physically reasonable singlet radiative rate of ~ 
10~® s~^ (see, for example, (49|). 

Fig. [S] shows the CIS and exact calculations of 
the radiative decay rates for the lowest singlet and 
triplet states. The CIS approximation does a rea- 
sonable job of reproducing the exact results, al- 
though the relative errors in this experimentally 
observable factor are a little larger than those in 
the calculated energies because fc^ depends approx- 
imately on the inverse of the fourth power of en- 
ergy differences^ Further, for — — the 
\^ M LCT^) character of the first triplet state, and 
therefore its radiative rate, is exactly zero. Hence, 
the hybridisation of the metal and ligand orbitals 
is a key determining factor of the phosphorescent 
properties of organometallic complexes. 

There is an order of magnitude change in the 
triplets radiative lifetime resulting from a 0.1 eV 
change in e, captured by both the exact and CIS 
solutions. This huge sensitivity of the radiative 
lifetime of the triplet is observed in phosphores- 
cent complexes after small chemical modifications 
Thus, the CIS approximation captures the same 
physics as the exact solution; the sensitive depen- 
dence of the radiative rate on e may explain the sen- 
sitivity of organometallic complexes to small chem- 
ical changes. 

4. Discussion 

We have shown above (see particularly Fig. [6]) 
that the model Hamiltonian investigated here pre- 
dicts that small changes in the parameters lead to 
large changes in the radiative lifetime. This cor- 
responds well with what it known experimentally, 
where single chemical substitutions are found to 
lead to dramatic changes in the photo luminescent 
quantum yield 01 • These changes are captured re- 
markably well by the CIS approximation. This sug- 
gests that the CIS approximation will be reliable 
and accurate as one moves to more complicated 
models, e.g., where more than one ligand plays an 
important role in the optical physics. 



iCf. Eqs. \TE\ and [19] and note that |Ti) is domi- 
nated by pLC^) (so Et^ ~ E^LC^) with a perturbative 
\^MLCT^) component, and that Esj^ji^Qrj^i ~ mlct^ > 
so i^MLCT^Ti) oc (i?Ti -i?iA/LCTi)"'- 
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Given the accuracy of the CIS approximation it 
is interesting to briefly compare and contrast the 
CIS approximation with those made in TDDFT cal- 
culations. CIS is quite analogous to (and proba- 
bly slightly worse than) the approximations typ- 
ically made in TDDFT calculations. CIS starts 
from a Hartree-Fock, i.e., mean field, description of 
the ground state. The local density approximation 
(LDA) is also a mean field theory, albeit a slightly 
different flavour. The generalized gradient approxi- 
mation (GGA) and hybrid functionals include some 
corrections to the LDA. Therefore although DFT 
is in principle exact, as currently practiced, it re- 
mains close to a mean field field theory. To find the 
excited states one could solve the time-dependent 
Kohn-Sham equations [H^. However, it is often 
more convenient to make a linear response approx- 
imation, such as the random phase approximation 
or the Tamm-Dankoff approximation (TDA). The 
TDA is exactly equivalent to a configuration in- 
teraction singles (CIS) calculation |4^]. Thus our 
CIS calculation using a mean-field ground state is 
close in both spirit and computational detail to a 
TDDFT calculation using the TDA. This means 
that a comparison of the CIS results from our model 
Hamiltonian with its exact solution provides some 
insight into, and perhaps and approximate lower 
bound on, the expected accuracy of TDDFT calcu- 
lations for organometallic complexes. 

5. Conclusions 

The results presented above show that, for the 
lowest excited states, there are no significant dif- 
ferences between the exact and approximate wave- 
functions and energies. This means that the key 
radiative properties of the complex are accurately 
reproduced in the CIS approximation. These re- 
sults also show that there are broad ranges of pa- 
rameter values for which these approximations can 
reproduce the energies of the higher excited states. 
The errors introduced by making the CIS approxi- 
mation mostly effect the high energy states, which 
do not play a significant role in, e.g., the phospho- 
rescence. Hence, we conclude that electronic corre- 
lations do not play an important role in determining 
the optoelectronic properties of these complexes. 
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Figure 1: [Colour Online] (a) Singlet energies as a function of 
e and (b) the relative error in those energies on a logarithmic 
scale. In panel (a), the thick solid lines are the exact results, 
the points are the results of the CIS calculations, and the 
thin black lines are the = t^ =0 energies. These figures 
show that the energies of the ground state and singly excited 
states are well approximated in this parameter regime. The 
doubly excited states energies are most accurate away from 
avoided crossings. The labels on the right hand side of (a) 
indicate the = = eigenstates. The energies of these 
states are all linearly dependent on e, with gradient (riM)- 
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Figure 2: [Colour Online] Exact and CIS triplet energies as 
a function of e, with the exact singlet energies show for com- 
parison. The thick solid lines are the exact results for the 
triplets, the points are the results of the CIS calculations for 
the triplet states, the thin black lines are the = = 
triplet energies, and the dashed lines are the exact singlet 
results. This figure shows that for e > 0.25 eV, the low- 
est excited triplet state is nearly degenerate with the lowest 
excited singlet, and both are predominantly metal-to-ligand 
charge transfer states; while for —0.75 < e < 0.25 eV the low- 
est triplet is mostly x. The third triplet is nearly degenerate 
with the fourth excited singlet as both are predominantly 
\MH^) (and the only singlet— triplet splitting in this model 
occurs due to mixing with the \LC^) states). The labels on 
the right hand side indicate the = = triplet eigen- 
states. The energies of these states are all linearly dependent 
on e, with gradient (tim)- 
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Figure 3: [Colour Online] The error in the mean field wave- 
functions, defined as 1 minus the overlap between the exact 
and CIS states, versus e. The ground state and first two ex- 
cited states have errors that are typically much less than 1% 
for the parameter values in this figure; the inset shows the 
error in these states on a logarithmic scale. The largest error 
is in the three doubly occupied states around e = —0.75 eV, 
coinciding with the avoided crossing between the first two ex- 
cited states (see Fig. [T]|. The two other local maxima in the 
error occur around e = eV, where the three = 2 states 
have an avoided crossing, and at e = 1.1 eV, where \^LC'^) 
and \^ M LCT'^) have an avoided crossing. The errors in the 
wavefunctions for the ground state and first excited state 
are below 1% for the parameter values in this figure. Away 
from e = 1.1 eV, the error in the wavefunction of the second 
excited state is also below 1%. From this data we conclude 
that the approximate ground state and singly excited states 
will typically do quite well, while the higher excited states 
only do well away from avoided crossings, that is, when they 
can once again be well represented by the = = 
eigenstates. 
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Figure 4: [Colour Online] The mag nitude of the = = 
eigenstates in the exact (solid lines) and CIS (markers) third 
excited singlet state, S3 (green triangles in Figs. [T]and[3]|. 
The |0) and \^ M LCT^) components are not shown as they 
are not significant in the region plotted. Around e = there 
is a three-way avoided crossing, involving all three = 2 
states: I^LC^}, \'^MLCT^), and \^MH^). At e = 1.1 eV 
this excited state changes character from mostly | ^ MLCT^ ) 
to mostly \^LC^). In the CIS solution these transitions are 
sharp crossovers rather than the gradual crossings seen in 
the exact solution. Around e = —0.75 eV one can see the 
effects of avoided crossing between the \S\) and [52) states 
(composed primarily of the \^LC^) and \^MLCT^) states). 
The CIS procedure optimizes the energy of the \S\) state; 
at the avoided crossing where this state is not well described 



by a single t 



■ eigenstate the CIS procedure causes 



the higher excited states to also be strongly mixed. 
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Figure 5: [Colour Online] (a) Singlet energies as a function of 
l/jv/, and (b) the error in the CIS wavefunctions on a logarith- 
mic scale, over a large range oIUm- In (a), the coloured solid 
lines are the exact results, and the points are the results of 
the CIS calculations. It is clear that for Um ^ 5.5 eV, when 
the ground state is no longer predominantly the reference 
basis state, the mean field approximation breaks down. For 
Um > 5.5 eV, the ground state contains a large \^ M LCT^) 
component, and therefore is no longer closed shell. In this 
regime electronic correlations become important and the CIS 
approximation fails dramatically. 
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Figure 6: [Colour Online] The radiative decay rate of the 
lowest triplet state, as calculated by perturbation theory: 
in the exact solution (dotted line), and in the mean field 
solutions (blue crosses). The solid line is the radiative decay 
rate of the exact solution of the first excited singlet state, and 
the red crosses are the mean field solution. Note that there is 
no qualitative difference between the triplet radiative decay 
rate found from the exact and mean field solutions. The 
singlet radiative rate calculated from the mean field state is 
within 5% of the exact rate in this figure. 
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